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@ Introduction 


> The multiple gravity-assist, low-thrust (MGALT) model 
> Time regularization 
> multiple gravity-assist, low-thrust Sundman (MGALTS) model 


@ Universal Anomaly Algorithm 


> The Sundman transformation 
» Extending MGALT to MGALTS 
> The Multiple Gravity Assist Low-Thrust (MGALT) model 


e@ True Anomaly Algorithm 


> The generalized Sundman transformation 
> Extending MGALT to MGALTS 
> The Multiple Gravity Assist Low-Thrust (MGALT) model 


@ Analytic Derivatives 
@ Application: Rendezvous with Comet 45P /Honda-Mrkos-PajduSakova 


e@ Summary 
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Integral aspect of preliminary mission design 


e@ Provide an accurate mass budget, if not trajectory 
e EMTG, MALTO, PaGMO, GALLOP 
> Sims-Flanagan transcription, (Sims and Flanagan, 1999) 
e@ Medium-fidelity solutions used as initial guesses for flight-qualified solutions 
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Forward propagation 


Keplerian propagation mw 
Mission clock time flow Backward propagation 
—> 


e@ Sims-Flanagan and patched conic approximation, phase divided into N 
control segments 


e@ Control nodes evenly spaced in time At = = 
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e@ Equal temporal spacing of control nodes is non-optimal for increasing 
eccentricity. 


e@ Apoapse raising is most efficiently achieved via maneuvering at periapse. 
e Using additional segments results in an increase in problem dimension. 
e Redistribution of the same number of nodes would be more desirable. 


@ This is typically achieved with a form of time regularization. 
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e Differential transformation introduced by Karl Sundman to regularize the 
three-body problem. 


e@ A new variable y replaces time as the independent propagation variable: 


dt = x (2) 


@ This results in an even geometric distribution of control nodes around the 
orbit. 
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@ The MGALT model encodes the phase flight time At, as a decision variable 


e@ MGALTS adds a total angular propagation variable Ayp 
e@ Each segment is propagated for a portion of the total independent variable: 
Axp 
Ax = 4 
e@ The propagation time associated with each Ay is computed using a universal 
Kepler solver in an inverse sense. 


> x is transcendental in Kepler’s equation, but is now selected by NLP solver 
==> iteration is no longer required in Kepler solver 


@ An additional nonlinear match point constraint is required to ensure match 
point epoch is continuous: 


rp—YrrF 

a= |e | 6 (3) 
MB— MNF 
tp —tpr 


e Analytic derivatives of the match point constraint are still possible to 
compute 


ee GALT Sun Ta 


e@ For MGALT, the maximum permissible Av in a segment is computed as 
follows: 


NactiveD TmaxAt 
Mk 


(4) 


Avmax = 


e@ For MGALTS, this can only be computed after the segment propagation time 
At has been determined 
@ Impulses shifted from the center of the segments to the boundary 


e@ The maneuver at the match point can potentially be twice the size of the 
others 


> This potential reduces as the number of segments increases. 
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@ In order to achieve a different control node distribution, consider the 
generalized Sundman transformation 


dt = c r7dx c€R and yEe{R:7>0} (5) 


@ Universal Kepler propagator based on {U;,(x; a@)}, derivations based on 7 = 1 
corresponding to the eccentric anomaly. 


e@ Require a relationship between Af and Ay 


2r,sin( SAf) 


VP cos( SA f)—oxsin( SAS) 5a = 0 
2 =] Var,sin( SA f) 
Ax = VYatan VP artad | mo 4 >0 (6) 
2 ae V—argsin( dA f) 
Jratanh E cos( SA f)—oxsin( FAS) 5 < 0 


Parameter Value 
Maximum allowed initial mass 4000 kg 
Launch window 5 years 
Maximum flight time unbounded 


Earliest allowed launch date 
Latest allowed arrival date 
Launch vehicle 

Maximum launch C3 

Launch declination bounds 
Thruster 

Throttle logic 

Thruster duty cycle 

Solar array BOL power (launch at 1 A.U.) 
Power margin (dpower) 
Post-launch checkout coast 
Number of segments per phase 
Ephemeris 

Objective function 


August 27°", 2016 

January 1°*, 2030 

Atlas V (555) NLS-2 
21.7156 km? /s” 

[—28.5°, 28.5°] 

2 x NEXT TT11 high thrust 
minimum number of thrusters 
0.95 

40.0 kW 

15% 

60 days 

200 

SPICE 

maximize final mass 


e Solar electric propulsion model adheres to Discovery class mission proposal 
requirements 
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Event # 2: 
Lt rendezvous 
7 11/21/2017 


45P , 
/ Cy = 21.716 km?/s? 


12/29/2027 
m = 3238 kg \\ DLA = 20.1° 


m = 4000 kg 


1e8 


Figure: Rendezvous with 45P using MGALT with 200 segments, Atmission = 3690 days 
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Event # 2: 
Lt rendezvous 
45P 
10/28/2027 
m = 3217 kg 


> [ Event # 1: 
Launch 

Earth 

11/10/2019 

Cy = 21.716 km2/s? 


DLA = 16.6° 
m = 4000 kg 


g 


384 


= 


Figure: Rendezvous with 45P using MGALTS with 200 segments, Atmission = 2909 days 
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— Distance from SUN (AU) 
-- Control magnitude 


@ MGALT requires an 
additional revolution of 
the sun 


SFE LO VEO BT VENT TOPO WET OEOTI WT RISOTTO @ Additional nodes at 

periapse affords 

Or ane MGALTS more control 
authority 


e MGALTS TOF is 
substantially less 


01-072020 1-0-2021 ——OFOT-2022—OLOT SOTO I0ad OTS 0-01-2036 OT-OTa0aT 
Epoch 
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Event # 2: 
Lt rendezvous 


5/19/2028 
m = 3138 kg 


2 (Event #1: 
Launch 


Earth 
11/16/2018 


= 


g 
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Figure: Rendezvous using MGALT with 400 segments, Atmission = 3472 days 


e@ MGALTS can achieve the same transfer using far fewer control nodes 
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Event # 2: : Event # 1: 
Lt rendezvous a = Launch 

45P < Earth 

10/29/2027 Z \ 11/8/2019 

m = 3217 kg Y \ Cy = 21.714 km?/s? 


DLA = 16.7° 
m = 4000 kg 


Figure: Rendezvous using the MGALTS true anomaly-based transcription with 200 
segments, Atmission = 2913 days 
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Event # 1: 
Launch 


an e@ Using true anomaly as 

Cy = 21.716 km?/s? . 

0] DLA = 13.3° the independent 
variable, nodes are now 


5 concentrated at 
a periapse 
/ Ss @ Short flight time 
pump-up maneuvers are 
accessible 


e Useful inside an outer 
loop that is capable of 
exploring wide flight 
time ranges 


1e8, 


Figure: MGALTS true anomaly, 200 segments, 
Atmission = 2350 days 
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e@ A medium-fidelity low-thrust model has been developed that solves the 
time-regularized Kepler problem exactly. 


e Straightforward extension of the MGALT model in use by many preliminary 
design tools 


e Kepler propagation is made iteration free at the expense of one additional 
decision variable and constraint per phase. 
e Especially useful for trajectories requiring large eccentricity increases 
> e.g. Pump-up/down maneuvers 


e@ MGALTS can perform certain transfers using fewer control nodes than 
MGALT. 


@ Certain solution families can only reasonably be accessed using MGALTS. 
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Thank You 


e@ Ricardo Restrepo and Ryan Russell for inspiring the iteration-free Kepler 
propagator. 


e@ Kyle Hughes and Jeremy Knittel (GSFC 595) for reviewing. 
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